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We show that the shear viscosity of A(^* theory is sensitive at next-to-leading order to soft physics, 
which gives rise to subleading corrections suppressed by only a half power of the coupling, rj — 
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for the pressure. 



The series appears to converge about as well (or badly) as the series 
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I. INTRODUCTION 

For some time there has been a growing interest in the 
theoretical prediction of long-time dynamics and equili- 
bration in hot quantum field theories. One driving mo- 
tivation has been understanding the early universe, par- 
ticularly the possibility of electroweak baryogenesis (for 
an older but still quite pertinent review see [1]). In this 
case one is interested in dynamics of the electroweak sec- 
tor, SU(2)-fHiggs with a couphng a^, ~ 1/30. Another 
motivation is attempting to describe high energy heavy 
ion collisions, which are being experimentally probed at 
RHIC and will soon be probed at even higher energies at 
the LHC. Since one of the original motivations for this 
program was to explore QCD in a regime where it is ther- 
mal and (relatively) weakly coupled, it makes sense to see 
how far one can go in its description using weak-coupling 
tools, even though the gauge coupling is not expected to 
be very small, as ^ 1/3 at best. 

Some of the most interesting and difficult problems 
involve understanding long-time dynamics and the ap- 
proach to equilibrium. Transport coefficients are a theo- 
retically well-defined and interesting subset of these prob- 
lems. Recently there has been great interest in the shear 
viscosity, because RHIC results on elliptic and radial flow 
and spectra seem to indicate that the viscosity is surpris- 
ingly small [2-4]. On the theoretical side, great strides 
have been made in the theoretical evaluation of transport 
coefficients. In a pioneering paper, Jeon showed how to 
calculate the shear viscosity of A^"* theory at leading or- 
der in the coupling [5]. He showed that the evaluation 
reduces to a problem in kinetic theory. Since then, ki- 
netic theory treatments of QED and QCD transport co- 
efficients (electrical conductivity, conserved number dif- 
fusion, and shear and bulk viscosity) have appeared at 
leading- log [6] and full leading order [7, 8]. The ki- 
netic treatment has also been justified diagrammatically 
at leading-log [9, 10] and leading order [11]. 

However, except for two unrealistic simplifying limits 
(large A^f gauge theory [12] and large N 0{N) scalar field 
theory [13]), there are no results beyond leading order 
in weak coupling for any long-time dynamical quanti- 
ties in any interacting relativistic 3-1-1 dimensional gauge 
theory [19]. We feel this is a big hole in the literature, be- 
cause without any knowledge of subleading corrections. 



it is hard to know how quickly perturbation theory con- 
verges. On the other hand, a subleading result gives at 
least a hint of how fast the series converges; when the 
subleading correction becomes as large as the leading or- 
der result, perturbation theory has almost surely failed. 
With this in mind, this brief report will compute the 
shear viscosity of the simplest toy theory, N component 
X(j)'^ theory at general N, at next-to-leading order. 

The reason this is possible, and a main reason it is in- 
teresting, is that the subleading correction emerges from 
soft physics. In this theory, the plasma of scalars in- 
duces an effective mass for all scalar excitations, m^j^ — 
^^XT"^. But soft scalars play an important role as the 
"targets" in scattering processes. The scalars which are 
significantly influenced by this effective mass represent 
the target in 0{mth/T) of scattering events. Therefore 
the scattering rate, and with it all long-timescale dynam- 
ical quantities, receive 0{VX) corrections. Actually, both 
this thermal mass and the mass dependence of the shear 
viscosity were mentioned by Jeon [5], so his 1995 paper 
technically represented an NLO calculation, not just the 
leading-order one he claimed. But this was not men- 
tioned (or recognized) in that paper, so we feel that a 
clear discussion is warranted. 



II. LEADING-ORDER REVIEW 

First we review Jeon's leading order calculation, em- 
phasizing where higher order corrections can arise. The 
shear viscosity is deflned as a correction to the (local) 
equilibrium form for the stress tensor when a fluid un- 
dergoes nonuniform flow. For a flow velocity Vi satisfying 
djVj — 0, the stress tensor in the local rest frame is 



diV 



+ Oid'), (1) 



with 77 the viscosity. A fluctuation-dissipation relation 
(Kubo formula) relates this coefficient to a certain stress- 
stress correlation function [15]: 
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(2) 



where TTy = Ty 



j-Tkk- This relation is the starting 
point for a perturbative analysis. One wants the correla- 
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FIG. 1: Two leading order diagrams contributing to shear 

viscosity. Blobs arc strcss-tcnsor insertions and heavy lines 
are self-energy resummed as shown. 



tion function of two stress operators, at zero spatial mo- 
mentum and in the limit of vanishing frequency. Figure 1 
shows the obvious leading-order diagram, together with 
another diagram which is also leading-order for a subtle 
reason. The heavy lines represent self-energy resummed 
propagators. The propagators with the same label carry 
opposite momenta for all values of all loop momentum 
integration variables; therefore when one goes on-shell, 
the other goes on-shell simultaneously, leading to a pinch 
singularity. If not for the width F, the result would di- 
verge; instead it is enhanced by a factor of 1/F ~ 
Therefore we must sum all diagrams with n such pinches 
and 2n — 1 loops. More loops without extra pinches are 
generally suppressed by a power of A per loop. 

The set of diagrams which we must resum turns out 
to be all diagrams structurally like the second diagram 
of Figure 1, and up to subleading in A corrections, their 
resummation gives precisely the same description of ij 
which we would get from kinetic theory. The kinetic de- 
scription takes the plasma to be made up of weakly cou- 
pl(;(l quasiparticles which undergo occasional scatterings 
described by a Boltzmann equation. The relevance of 
quasiparticles arises because of the importance of nearly 
on-shell propagators in evaluating the diagrams; the scat- 
tering processes arise from evaluating the self-energy cor- 
rections to these propagators, which give rise to F. Ac- 
cording to kinetic theory, the stress tensor is determined 
by the distribution function for particles /(a;,p) through 



Tij{x) 



d^P PtPj n, ^ 



(3) 



This distribution function evolves according to the Boltz- 
mann equation. 



^/(a:,p) Pi dfix,p) _ 

dt dxi 
1 f _cfkd^p'd^k 



(4) 
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x(/(p)/(k)[l+/(p')][l+/(k')] 

-/(p')/(k')[l+/(p)][l+/(k)] 



In writing Bose stimulation functions we have assumed 
all bosonic species; wc have also neglected higher order 
scattering processes, as justified by our expansion in A. 

Nonuniform flow velocity means that the equilibrium 
distribution function would be 

/o(a;,p) = (exp[(pO - ViPi)/{Ty/l^)] - l)"' 

and the spatial derivative in Eq. (4) acts on this to give 
{PiPjdiVj)fo{l+fo)/{p°T). This term drives the system 
from equilibrium. The collision term vanishes if / = /o 
but is linear in departures Sf ^ f — fo] such a depar- 
ture will arise and grow until the collision term and the 
gradient term cancel (remember that we want the van- 
ishing frequency or long-time behavior). The departure 
from equilibrium will be of the same structural form as 
the spatial derivative term; 



f{x,p) = fo 



idiVj)^Ml+fo)x{\p\), 



(6) 



with x(|p|) a function of the magnitude of p only, de- 
termined by solving Eqs. (4,5). This is most con- 
veniently done by defining an inner product {f\g) = 



/ ■(2^fiP)9{p) and a functional Q[x] [6], 

Q[x]=T{x\^-^}-hx\Cnn\x) 



(7) 



where Ciin is a linearization of the collision operator given 
below. The x which extremizes Q solves the Boltzmann 
equation, and the extremal value of Q is 15r]/2. 



III. SUBLEADING CORRECTIONS 

The advantage of the diagrams is that they help us 
see where the perturbative argument that further correc- 
tions are down by 0(A) can break down. Specifically, for 
generic momenta an extra loop is 0(A). The self-energy 
resummation is an 0(A) correction imless the propaga- 
tor momentum is soft or close to the light-cone. Near the 
light-cone, the width correction is essential, but it is al- 
ready incorporated in the leading-order calculation and it 
naively receives only 0(A) corrections from higher order 
diagrams. The one-loop self-energy leads to a mass cor- 
rection mentioned in the introduction, m^j^ = ^^^j^AT^ 
for N component scalar theory with 0{N) symmetry and 
interaction ■^{4>n4'n)^- This mass correction shifts the 
propagating pole by 0(A) for generic momentum p ^ T. 
However, for p ~ rrith ~ TV NX, the shift to the propa- 
gating pole is an 0(1) correction. Therefore, in any di- 
agram where one of the lines is soft, p ^ mth, there will 
be 0(1) interaction corrections not taken into account in 
the leading-order analysis. 

How important are soft momenta? It turns out that 
x{p) vanishes as x{p) ^ P at small momenta, meaning 
that particles at these momenta are very near equilibrium 
(5) and do not contribute significantly to Tij — P6ij , despite 
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the /o(l+/o) enhancement factor. Therefore the first 
term in Eq. (7) is only sensitive to dispersion corrections 
at 0(A). However, such particles are of some importance 
as scatterers. To see this, first let us estimate what frac- 
tion of all particles are so soft. The total particle number 
density is ntot = N J{d^p/{2Tr f)fo - NT^. The soft 
contribution is ngoft ~ N J {d'^p/ {27T)^)foQ{mth—p) ~ 
NttlIy^T ^ (7VA)ntot- (Phase space gives three factors 
of mth but the distribution function is /o ^ T /mt]^ in 
the soft region.) Soft particles are therefore an 0(NX) 
fraction of all particles. In addition, scattering cross- 
sections gcnerically scale as cr 1/s (the Mandelstam 
variable), and when a. k = 0{T) particle scatters from 
a p = O(mth) particle, s ~ —2pnk^ ~ mthT is smaller 
than s T"^ valid for generic scattering. Therefore the 
cross-section is enhanced by a T/rrith factor, and the rel- 
ative rate of sca tteri ngs from a soft target is actually 
0(mth/T) = 0(VNX). 

So soft scatterings are an 0{VNX) fraction of all scat- 
terings. Do they matter to equilibration? The answer 
is yes; while the incoming particle is soft, most of the 
phase space for the scattering shares the final state ener- 
gies roughly equally. This means that the hard particle 
retains its direction of flight but substantially changes its 
energy, and since the departure from equilibrium has a 
strong energy dependence, this is relevant. To calculate 
it, we need to perform the scattering calculation more 
carefully, incorporating the soft phase-space dispersion 
relation. 

The collision term we need for the functional Q[x] is 



ixlCM^Jj^ 



d^pd^kd^p' d^k' 
)i2 2p0 2fc0 2p'0 2fc'0 T5 

x{2n)^d^{p+k-p'-k') 
1, 



(8) 



Xijip) 



x^\Mffoip)foik)[l+foip'm+Mk')] 

X iXij{p)+Xi]ik)-Xij{p') - Xtjik')f , 

PiPj - ^ijP^ 



xip)- 



(9) 



The phase space is most conveniently dealt with in terms 
of the total energy iv = p°+k^ and momentum q = |p+k|, 
two particle momenta, and an azimuthal angle. For mass- 
less dispersion relations, this "s-channel" parametriza- 
tion reduces to [12] 



d^pd^kd.^p' d^k' 4 4 , , 

— -o(27r) 6 {p+k-p-k 



(27r)i2 2p" 2A;0 2p'" 2/e' 



2tt 



(10) 



The relations between these variables and the relative 
angles of the particles are given in [12], and can be used 
to evaluate Eq. (8). 

The particle energies are p, p', co—p, and cu—p'. The 

u! ^ m region is to** suppressed, as all four variables 
must be ~ m. [There is a (T/m)^ enhancement from the 



population functions, but a {m/T)'^ suppression from the 
small value of Xij^ so this region really is 0{X^).] A re- 
gion with one soft particle is suppressed; for instance, 
for soft p, both p and uj—q must be 0{m). However, one 
of the statistical functions in Eq. (8) is 0{T/m) in this 
case, and there is no suppression from the values of Xij- 
Therefore the contribution is 0{m/T), as stated above. 
Note however that two soft particles (say, p,p' ~ m) is 
0{'m^) suppressed, compensated by two powers of popu- 
lation functions but further suppressed by cancellations 
in the % functions; so this region is 0(A) at most. 

For the massive case it is more convenient to work 
in terms of particle energies rather than momenta. As 
just discussed, we need only consider the case where one 
particle is soft. In the small region, and taking w ^ 
mth, the p", q phase space is modified by 

dq dp= dp dq ^ dp° dq, (11) 

J Jq Jui-2p Jm Jui-p°-p 



with p = \J (p^)"^ — m^jj. If we consider the difference 
between the collision integral with and without massive 
dispersion, we want the difference between these integra- 
tions. In this difference region, we can take p° small, 
allowing the substitution fo{p^) = ^ and f{k) = f{co). 
Furthermore, we may neglect Xijip) ^^'^ set Xij{k) = 
{kikj /uj)x{ijj)- In this region, k and the final state mo- 
menta p' , k' are collinear, so 

iXtj ip)+Xij (k) -Xtj {p')-Xi] {k')f 
- {^Xi^) -p'xip ) - {'^-p')x{'^-p')f ■ (12) 

The q,p^ and u^p' integrals then factorize, such that 



(x|c lin,m=0 

p pUJ rp f- f-UJ—%P-^'p rp 

dp" dq-- dp'> dq- 

Jo Ju-2p0 P Jm Ju>-p°-p P 

• I I ^ X ^ J^dp'Mu,)e^MpVo{^-p' 



8 287,-. JO 
X {ijxitj) - p'xip') - {LJ-p')x{iJ-p')f . (13) 

Here the factor of 4 counts the fact that there are 4 such 

"corners" of phase space where a particle becomes soft, 
and we have only considered the contribution of one of 
them. Note that the p' integral is insensitive to the edges 
p' ~ and p' ~ w of its range, so there is no problem of 
overcounting. The bracketed integral gives nmthT. 

We want to perturb in the small mass. Therefore it is 
fair to write the collision operator as 



Clin — Clin,m=0 + 



(14) 



with 5C equal to minus the expression in Eq. (13). The 
correction to Q in Eq. (7) is therefore — ^(x|^C|x). The 
change in the extremal value at linear order in 6C is just 
found by evaluating this using the unperturbed extremal 
value of X- 
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What remains is to solve the variational problem with 
the unperturbed collision operator and to substitute the 
result for x iiito Eq. (13). We handle the variational 
problem using the multi-parameter Ansatz method de- 
veloped in [6]. Since we are tuning the function % using 
the m = collision operator, the evaluation of (xl<^C|x) 
only improves linearly with the quality of the variational 
Ansatz. Therefore we must use a large number of the 
variational functions presented there; going from 1 to 2 
basis functions changes our result by 4.6%, going from 2 
to 4 functions changes it 2.3%, going from 4 to 6 func- 
tions changes it 0.13%, and going from 6 to 8 changes it 
0.01%. Our result in the one-component theory, using an 
8 basis function expansion, is 

_ 3033.54+ 1548.3 TOth/r ^3 

V — ^ 

where both integration and basis size errors should be 
smaller than the last digit shown. The viscosity increases 
because the thermal mass reduces the number of soft tar- 
gets, and viscosity is inversely related to the scattering 
rate. Jeon found a leading result 5% smaller and a mass 
correction 15% larger, presumably because of lower qual- 
ity numerics. The result can be extended to the 0{N) 
symmetric theory by scaling the number of species by 
N and the scattering rate and thermal mass each by 
{N+2)/3. 

IV. DISCUSSION 

We have shown that the first corrections to the leading- 
order calculation of the shear viscosity in scalar 0^ theory 



arise at order 0(\/A). They originate because scattering 
processes are sensitive to soft {p ~ T^/X) particles at 
that order, and soft particles' behavior is modified by 
plasma effects at the 0(1) level. All of these remarks 
should apply to gauge theories including QCD. However, 
the details of the calculation in those cases will be more 
complex, since the dispersion relations are modified in 
a more complex manner in gauge theories, and in non- 
abelian gauge theories the mutual interactions between 
soft particles are also strongly modified [16]. 

One of the motivations for this work was to analyze 
the convergence of the perturbative treatment. What we 
have shown is that the first subleading correction is a 
small correction provided mth ^ 2T. This is satisfied if 
A <^ j^q^ . The thermal pressure converges in a compara- 
ble sized range [17]. So does vacimm perturbation theory; 
the ratio of the /3 fimction to the coupling ^ = |g^A 

is small for A <C Therefore the subleading correc- 

tions only become large when the theory is approaching 
its Landau pole. Since the origin of the leading-order 
correction is well understood, we might hope that the ex- 
plicit inclusion of the thermal mass in the external state 
dispersion, and a resummation of the thermal mass via 
a gap equation, might improve the convergence of the 
series, as has been argued for the calculation of the pres- 
sure [18]. Aarts and Martinez- Resco have done this in 
the large N theory, where it is well justified at leading 
order in a 1/A'' expansion [13]. 
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